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Comprehensive Monte Carlo simulations of the short-time dynamic behaviour are reported for 
the three-dimensional Ising model at criticality. Besides the exponent 6 of the critical initial increase 
and the dynamic exponent z, the static critical exponents u and /3 as well as the critical temperature 
are determined from the power-law scaling behaviour of observables at the beginning of the time 
evolution. States of very high temperature as well as of zero temperature are used as initial states 
for the simulations. 
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I. INTRODUCTION 

Critical properties of many magnetic materials can be 
described by a simple Ising model, 



H = K S,Sj 



(1) 



where Si = ±1 represents the spin of site i, and the sum 
extends over nearest neigbours only. The factor l/fc^T is 
included in the coupling constant K. The Ising model has 
been solved exactly in one and two dimensions. However, 
for higher dimensions, there exist extensive perturbative 
analyses based on renormalization group methods and 
numerical investigations with Monte Carlo methods. 

It was traditionally believed that universal scaling be- 
haviour exists only in or near thermodynamic equilib- 
rium. Recently, it has been argued theoretically [Q that 
some dynamic systems already exhibit universal scaling 
behaviour in the macroscopic short-time region of their 
dynamic evolution. The main point is that universality 
and scaling emerge after a microscopic time scale tmic, 
which is sufficiently large in the microscopic sense but 
still short in the macroscopic sense. This statement has 
been proven valid in a series of numerical simulations of 
various statistical systems. More interestingly, as well as 
the new critical exponent 6, describing the critical initial 
increase of the magnetization, and the dynamic critical 
exponent z, short-time critical dynamics provides a mea- 
sure of the static critical exponents and of the critical 
temperature. This dynamic approach is free of critical 
slowing down since the spatial correlation length is still 
small within the short-time regime, even at or near the 
critical point^. 

Up to now, systematic numerical simulations of the 
short-time critical dynamics have been carried out 



mainly in two-dimensional systems 0. For the three- 
dimensional Ising model, the power-law decay of the mag- 
netization starting from an ordered state has been simu- 
lated [^j4j , the new exponent 6 and the dynamic exponent 
z have also been obtained However, a complete un- 

derstanding of the short-time dynamic behaviour of the 
three-dimensional Ising model is still necessary, since it 
is a very important model. Especially, a systematic test 
of the short-time dynamic approach to the determination 
of all the critical exponents and the critical temperature 
in a three-dimensional system is important. 

In this paper we report a comprehensive investiga- 
tion of the short-time critical dynamics of the three- 
dimensional Ising model. A power-law behaviour of the 
autocorrelation, the second moment and the Binder cu- 
mulant is observed, in addition to that of the magneti- 
zation. The results support fully the short-time dynamic 
scaling. For the first time we extract the critical tempera- 
ture and the static exponents v and (3 from the short-time 
scaling behaviour. Our results for the static exponents 
and critical temperature agree well with those obtained 
in the extensive studies at thermodynamic equilibrium. 

In the next section, a scaling analysis of the short-time 
critical dynamics is given. Numerical results are pre- 
sented in Sec. III. The last section contains the summary 
and discussion. 



II. SCALING RELATIONS 

Using renormalization group methods, Janssen, 
Schaub and Schmittmann have shown that far from 
equilibrium, in a macroscopic short-time regime of the 
dynamical evolution, there already emerges universal 
scaling behaviour in the 0{N) vector model. The re- 
laxation process considered is one of a system initially 
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in a high-temperature state, suddenly quenched to the 
critical temperature Tc and evolving with dynamics of 
model A. For an initial state with a non-zero but small 
magnetization (mg <^ 1), a generalized dynamic scaling 
form has been derived with an e-expansion for the 0{N) 
vector model, 



Thus it can be deduced that M^'^\t,L) ^ L""^, where d 
is the dimension of the system. Taking this into account, 
Eq. M) yields a power-law behaviour 



M^^Ht)^t'\ C2 - (d-'^^j -■ (8) 
M^'^) (t, T, L, Too) = 6-'='^/'' M^'') (b-H, b^/^r, b^^L, b^'^mo). An analysis of the autocorrelation (for toq = 0) 



(2) 

In Eq. (0), AI^''^ is the fcth moment of the magnetization. 



(3) 



t is the time of the dynamical relaxation, L is the lattice 
size. 



T-T, 
T, 



(4) 



is the reduced temperature, and 6 is a spatial rescaling 
factor. The quantity is a new independent critical 
exponent, the scaling dimension of the initial magnetiza- 
tion mo . The interesting and important point is that the 
critical exponents /3, v and z in Eq. (^) are exactly those 
usually defined in equilibrium. These exponents can thus 
be extracted from the short-time critical dynamics. 

Equation (||) provides the time evolution of the mag- 
netization, with k = 1. Taking b — t^/^ for the spatial 
scaling factor gives 

M(t, r, mo) = t^^/'^^ A/(l, t^/^^V, t^o/^mo) (5) 

The expansion has been performed with respect to the 
small quantity t^°/^mo and M(t, mo = 0) = has been 
used. It has been implicitly assumed that L is sufficiently 
large. Exactly at the critical point (r = 0), Eq. (H) pre- 
dicts a power-law behaviour of the magnetization in the 
short-time region. 



Xq 



(3\ 1 



(6) 



Up to now, analytical calculations for the 0{N) vector 
model and numerical simulations for a variety of statis- 
tical systems show that 9 > 0, i.e. the magnetization 
undergoes an initial increase. This is a very prominent 
phenomenon in the short-time critical dynamics. 

Now we consider the case toq = 0. Using Eq. (||) for 
the second moment of the magnetization at the critical 
temperature gives 



(7) 



In the beginning of the time evolution the spatial cor- 
relation length is still small, even at the critical point. 



A{t) 



N 



shows that it obeys a power-law 



A{t) - r 



d 

Ca = I 

Z 



(9) 



(10) 



In summary, simulations of the dynamic system start- 
ing from small or zero initial magnetization, at or near 
the critical point, allow the quantities 0, ci and to be 
measured and thus a determination of the critical expo- 
nents 0, (ijv and z separately. In principle, the critical 
temperature itself can also be determined from the loca- 
tion of the optimal power-law behaviour of the magneti- 
zation within the critical region |^,^ . However, a similar 
but more accurate determination of the critical temper- 
ature will be presented below. This is also the case for 
the determination of the exponent v from the derivative 
of M with respect to r. 

In the above considerations the dynamic relaxation 
process was assumed to start from a disordered state with 
vanishing or small magnetization mo. Another interest- 
ing and important process is the dynamic relaxation from 
a completely ordered state. The initial magnetization be- 
ing exactly at its fixed point mo = 1, a scaling form 

M^^) (t, r, L) = b-^^'" mC') {b-H, b^/^T, b-^L) (11) 

is expected. This scaling form looks the same as the dy- 
namic scaling form in the long-time regime, however, it is 
now assumed already valid in the macroscopic short-time 
regime. 

For the magnetization itself, b = i^/^ yields 



This leads to the power-law behaviour 



M{t) - t-'^S ci = ^ 
vz 



(12) 



(13) 



at the critical point (r — 0). For small but nonzero 
T, the power-law behaviour of the magnetization will be 
modified by the scaling function M(l, i^/'^^r), thus al- 
lowing for a determination of the critical temperature . 
Taking the derivative with respect to r on both sides of 
Eq. (|l|) and fixing b = t^l"" again, gives the logarithmic 
derivative of the magnetization 
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drln M{t,T)\r=o - t-''\ cn = —. (14) 

vz 

Here, unlike the relaxation from a disordered state, 
the average magnetization is not zero. A Binder cumu- 
lant V{£) can be obtained using the magnetization and 
its second moment. Finite size scaling shows that 

Thus, the short-time behaviour of the dynamic relax- 
ation starting from a completely ordered state is also suf- 
ficient to determine all the critical exponents /3, v and z 
as well as the critical temperature. In practical simula- 
tions, these measurements of the critical exponents and 
critical temperature are usually better in quality than 
those from a relaxation process starting from a disor- 
dered state. 



III. NUMERICAL RESULTS 

We have performed simulations on three-dimensional 
lattices of linear sizes L = 32, 64 and 128 (in a particular 
case L = 256), starting either from an ordered state or 
from a high-temperature state with zero or small initial 
magnetization. In the latter cases the inital magnetiza- 
tion has been prepared by flipping in an ordered state 
a definite number of spins at randomly chosen sites in 
order to get the desired small value of mg. Starting from 
the initial state, the system has been updated by a heat 
bath Monte Carlo algorithm. A unit in time is defined 
as a complete update of all the spins in the lattice. Sim- 
ulations have been performed, depending on the inital 
magnetization, up to i = 1000. The preparation of the 
inital magnetization and the update up to the maximal 
time has been repeated 1000 or 5000 times, depending 
on the lattice size (much less for L = 256), and the mag- 
netization, the second moment of the magnetization (|^) 
or the autocorrelation (^ has been measured. In gen- 
eral, several runs of this kind have been performed to 
estimate the statistical error. The simulations used the 
value Kc — 0.22166 for the critical point, taken from 
Ref. [|l^. Simulations have also been carried out in the 
neighbourhood of Kc to extract the critical point and the 
critical exponent v. 



valid only in the limit mo = 0. However, practical mea- 
surements can only be carried out for finite mo and the 
measured exponent 9 may show a weak dependence on 
mo- The data must thus be extrapolated for mo — 0. 

We have measured the magnetization M{t) for mo = 
0.02, 0.04 and 0.06 for the largest lattice, L = 128. For 
mo = 0.02 we have performed four runs, in the other 
cases three runs up to < = 300. Fig. |l| shows on a log-log 
scale a pronounced power-law behaviour starting at the 
onset of the time evolution. A finite mo effect is observed 
at the larger times. 
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FIG. 1. Time evolution of the magnetization M(t) for 
three values of the initial magnetization mo for L — 128. 

In order to determine more clearly the region where 
the power-law is strictly fulfilled, we divided the time 
scale into non-overlapping intervals Aq{t) = [t, 9^], with 
g < 2, and measured the exponent 6{t) separately in each 
interval. Since the intervals contain more points at later 
times where the fluctuations are larger, the errors in each 
interval are comparable. Already from the second bin the 
slope shows a stable behaviour. For mo = 0.02 slope is 
quite stable up to f = 100. Therefore 9 was calculated 
from a least-squares fit in the time interval [2, 100]. 



mo 





0.02 


0.04 


0.06 


6 


0.108(2) 


0.1059(20) 


0.1035(4) 


0.1014(5) 



TABLE I. The exponent 8 measured for L = 128 for 
different values of the initial magnetization mo. The value 
^(mo — 0) is the result of an extrapolation (see text). 



A. Evolution from a disordered state 



Magnetization 

At the beginning of the time evolution, for sufficiently 
small t^f'/^mo, the magnetization undergoes a power-law 
initial increase. For mo and t not too small, the power- 
law behaviour will be modified. This may be called 'finite 
mo effect'. The possibility of a finite size effect also has 
to be evaluated. Furthermore, scaling form (||) is strictly 
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Similarly, we chose the interval [2, 50] for toq = 0.04. For 
mo = 0.06 the restrictive interval [2, 15] was used. Ta- 
ble I presents the results. From these results, an extrap- 
olation to Too = yields — 0.108(2). Careful analysis 
of the data for L — 6A shows that the finite size effect in 
L = 128 is already negligible. 

Autocorrelation 

Measurements of the autocorrelation function suffer from 
large fluctuations in the region of large t, since it decays 
by nearly three orders of magnitude between t = and 
t = 100. Fig. ||(a) shows results for four runs for L = 128 
on a log-log scale. For the lattice size L = 64 three runs 
were performed. 




-1.2 - 



-1.6 - 
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FIG. 2. (a) Autocorrelation A{t) and (b) the slope Ca{t) 
of the autocorrelation function for L = 128, in the intervals 
Ai.sW = [i, 1.5*]. 

As described above, we have divided the time scale into 
non-overlapping intervals Ai.5(t) and calculated CaU) in 
each interval. This is exemplified for L = 128 in Fig. @(b). 
On this plot Ca{t) is rather constant after the first few 
bins, while large fluctuations start a.t t > 70. Thus, the 
interval t = [10, 70] is used for L — 128 and 64 in order 
to obtain Ca- The results are given in table With 
no apparent dependence on lattice size, the two results 
were simply averaged to estimate Ca for L = oo (second 
column). 

Second moment of the magnetization 

The second moment of the magnetization has been sim- 
ulated for L = 128 in three runs up to t = 512, and for 
L = 64 in three runs each up to t = 100. Fig. ^ shows 
M(2)(i) for L = 128 on a log-log scale. 
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FIG. 3. Time evolution of the second moment of the mag- 
netization, M'^'(t), for L = 128. 

This curve as well as the slopes in bins show slightly 
larger fluctuations at higher i, but no departure from a 
power-law. Therefore, the full intervals t = [2, 512] for 
L = 128 and t ~ [2, 100] for L = 64 were used to com- 
pute the values of C2 reported in table ||. For L = oo, 
the estimate C2 = 0.970(11) provided is obtained by av- 
eraging the results for L — 64 and L — 128, since no 
dependence on the lattice size can be seen. 



B. Evolution from an ordered state 

For the determination of the critical exponents and 
critical temperature, the dynamic relaxation process 
starting from an ordered state has been proven advan- 
tageous over that from a high-temperature initial state. 
Indeed, the magnetization only decreases slowly in time 
from Too = 1. Furthermore, the magnitude of the magne- 
tization in the short-time regime is large, therefore statis- 
tical fluctuations are less prominent. We have measured 
the magnetization M(t) and the second moment M^'^\t) 
for lattice sizes 32, 64 and 128, and, with less statistics, 
for L = 256. Measurements have been performed up to 
t ~ 512 for L = 128, and up to i = 1000 for the magne- 
tization only. For the large lattice size L — 256 an av- 
erage over only 20 samples has been performed whereas 
for L = 128 and 64, over 1000 samples were taken for 
each run, and two to four runs performed in order to es- 
timate the statistical error. For L = 32, 5000 samples 
were collected. 

Magnetization 

Fig. U shows the evolution of the magnetization for all 
four lattice sizes on a log-log scale. It is interesting that 
for L = 256 down to L = 64 the measurements com- 
pletely overlap up to t = 1000. From Eq. (p^), the slope 
of the curves provides a measure of (3/vz. However, care- 
ful analysis reveals that the slope decreases weakly with 
increasing time. This suggests that a correction to scaling 
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should be considered in order to obtain accurate results. 
In comparison to the two-dimensional Ising model, the 
correction to scaling seems somewhat larger in three di- 
mensions [||. We have found that the best correction to 
scaling is given by the Ansatz 



M{t) = a r"^ e"^* ' . (16) 



Least-squares fits were performed in the interval t = 
[1,1000] for L = 64, 128 and 256. The resulting fit pa- 
rameters, averaged over the three lattice sizes and over 
the two runs for each lattice size, are ln(a) = 0.0033(10), 
ci = 0.2533(7), 7 = 0.1874(19), and 5 ^ 0.479(39). Ta- 
ble |l| reports the values for ci. Calculated values for 
M{t) using ( p^ ) are also provided on Fig. ^ (points) and 
show good agreement. Without the correction to scaling, 
the measured exponent ci would be one to two percent 
smaller. 
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FIG. 4. Time evolution of the magnetization M{t) for 
mo = 1 for L = 256, 128, 64, and 32. The points are obtained 
from the best fit to (^). 



Binder cumulant 

The second moment of the magnetization has been mea- 
sured up to i = 512 for L = 128, and up to t = 1000 for 
L = 64 and 32. A plot of the Binder cumulant Uit), de- 
fined in Eq. (p^), is shown on log-log scale in Fig. ^. The 
curves show that for L = 32 the power-law behaviour 
prevails only up to t ~ 100, while for the larger lattices 
it remains up to the maximal time. 
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FIG. 5. Time evolution of the Binder cumulant U{t) for 
mo = 1 and lattice sizes L = 128, 64, and 32. 

An analysis of the slope measured in the intervals 
^i.5(i) shows that the exponent cu{t) can be obtained in 
the interval i = 30 up to the maximal time for L = 64 and 
128. For L = 32 there is a plateau only over t — [10, 100]. 
Hence the results for this lattice size are not reported. 
The slope cu is calculated in the intervals t = [30, 1000] 
for i = 64 and i = [30, 512] for L = 128 and given in 
table ^ Within errors they are consistent and their aver- 
age is reported for L = oo in table ||. The error estimates 
are only based on statistics over a limited number of runs 
and exclude possible systematic contributions. 
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FIG. 6. Time evolution of the magnetization M{t) 
for mo = 1 for three values of the inverse temperature 
Ki = 0.22066, K = 0.22166, and K2 = 0.22266. 

Critical temperature and Cfi 

Figure ^ shows the magnetization M{t) for mo = 1 for 
three values for the inverse temperature Ki = 0.22066, 
K = 0.22166 and K2 = 0.22266, on a log-log scale. 
The critical point Kc can be estimated by searching 
for the best power-law behaviour with K between Ki 
and K2. Namely, the best straight-line fit to curves ob- 
tained by quadratic interpolation for Ki < K < K2 is 
sought. Restricting the interval iot = [30, 512], the result 
Kc = 0.22170(4) is obtained. The quadratic interpola- 
tion also provides the logarithmic derivative of the mag- 
netization with respect to t in Eq. (p^. This is shown 
in Fig. 1^. The slope over the interval [30, 300] provides 
C£i — 0.774(1). The errors are based on the statistics of 
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comparing only two data sets. A correction to scaling has 
not been considered. This could result in slightly larger 
errors. 



10 



3,ln M(t) 




FIG. 7. Logarithmic derivative of the magnetization with 
respect to r, obtained from a quadratic interpolation between 
the three curves shown in Fig. ^ taken a,t K = 0.22166. 





L — oo 


L = 128 


L = 64 


Ca 


1.362(19) 


1.358(19) 


1.366(20) 


C2 


0.970(11) 


0.966(16) 


0.975(3) 


Cl 


0.2533(7) 


0.2539(1) 


0.2527(8) 


cu 


1.462(12) 


1.453(17) 


1.471(4) 




0.22170(4) 


0.22170(4) 




cn 


0.774(1) 


0.774(1) 





TABLE II. Exponents Ca of the autocorrelation function 
and C2 of the second moment of the magnetization for dif- 
ferent lattice sizes obtained from a disordered initial state 
for mo — (upper part). Exponents ci of the power- law 
decrease of the magnetization and cu for the cumulant, as 
well as inverse temperature Kc and exponent cn for the log- 
arithmic derivative of the magnetization measured from an 
ordered initial state (lower part). For ci, the additional re- 
sult Cl = 0.2534(4) obtained from two runs with L = 256 has 
been included (see text). 



C. Critical exponents 

We now proceed to determine the critical exponents. 
Using the relaxation from a high-temperature state, we 
have measured the exponent 6 directly from the critical 
initial increase in magnetization and the exponents 
C2 and Ca from the power-law behaviour of the second mo- 
ment (H) and the autocorrelation (^^. Similarly, from 
measurements of the power-law behaviour with an or- 
dered initial state we have obtained ci, qi, and cu from 
the magnetization (|l3|), its logarithmic derivative (|l^ ) 
and the Binder cumulant (15). At this point, it is useful 
to recall the scaling relations: 



C2 



Cl 



Cu 



Cll 



A. 

vz ' 



_l_ 

vz ■ 



Our measurements of the exponent c\ are very accu- 
rate, while actual errors on and qi could be somewhat 
larger than those given in table ||. The value obtained 
for Cl agrees well with previous measurements [^IJ], al- 
though the approaches to correct for scaling differ. Our 
procedure is essentially similar to that used in Ref. [ pd] |. 
The dynamic exponent z can be estimated independently 
from C[/. Since the Binder cumulant ( p^ ) is constructed 
from the magnetization and its second moment, the es- 
timates are usually better than for a Binder cumulant 
constructed from the second and the fourth moments, 
typically used in the relaxation from a disordered state 
or in equilibrium. The exponent z can also be extracted 
from Ca + (^ and C2 + 2ci. Table III lists the results. These 



concur within statistical errors. Table IV reports the av- 
erage of the three values. With z, and from ci and qi, 
the exponents ^jv and v can be estimated. Results for 
all critical exponents and the critical point are given 
in table iV, 





djz 


z 


Ca +6* 


1.470(13) 


2.041(18) 




1.462(12) 


2.052(17) 


Cl + 2ci 


1.4766(8) 


2.032(11) 



TABLE III. Determination of z from three independent 
measurements of djz. 
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Our final result for the exponent given in table |, 
9 = 0.108(2) is consistent with an early estimate using a 
small lattice g but slightly larger than 9 = 0.104(3) ob- 
tained from damage spreading [|6| . It is not clear whether 
this difference comes from statistical or systematic er- 
rors. From equilibrium dynamics, the dynamic expo- 
nent z — 2.04(3) has been extracted ||l^. Values of 
z = 2.032(4) § and z = 2.04(1) [§ have been obtained 
recently from damage spreading and of z = 2.05(2) ||] 
and z = 2.04(2) [Q from large scale simulations of criti- 
cal relaxation, starting from an ordered state. Our value 
z = 2.042(6) is consistent with all of them. For the ratio 
(j/v and V, recent analytic calculations on the base of a 
renormalization-group expansion (in equilibrium) yield 
P = 0.327, P/v = 0.518 and ly = 0.631 |ll-|li, while nu- 
merical investigations in equilibrium by Ferrenberg and 
Landau yield /3 = 0.3258(44), = 0.518(7) and 
ly = 0.6289(8). Ref. [|oj provides a good review of earlier 
numerical values. Investigations using the cluster algo- 
rithm have yielded ly = 0.6301(8) and /3 = 0.3267(10) 0] 
and f3 = 0.3269(6) 0. In 111, the value = 0.6250(25) 
has been extracted from measuring the interface energy. 
Our value oi = 0.517(2) is accurate but the value of 
V is slightly larger. Our estimate of Kc is consistent with 



Kc = 0.2216595(26) in Ref 
ll, or = 0.2216544(3 
been quoted in Refs. pSl pOl 




Kc = 0.2216546(10) 
Similar results have 



e = 0.108(2) 



z = 2.042(6) 



(5/v = 0.517(2) 



u = 0.6327(20) 



(3 = 0.3273(17) 



Kc = 0.22170(4) 



0.104(3) 



12), 2.032(4) §, 



2.04(3) 

2.04(1) 1^2.05(2) 
2.04(2) 



0.518(7) liqi, 0.5185(16) |1 



0.6289(8) |1C , 0.6250(25) [JJ 

0.6301(8) M_ 

0.3258(44) Jl^, 0.3267(10) |l| 
0.3269(6) 



0.2216595(26) 
0.2216546(10) 
0.2216544(3) 



lei 



TABLE IV. Final results for all critical exponents and the 
critical point K^ (left). Results from earlier investigations 
discussed in the text have been collected in the right column. 



IV. SUMMARY AND DISCUSSION 

In the previous sections we have reported comprehen- 
sive Monte Carlo simulations of the short-time critical 
dynamics for the three-dimensional Ising model. Start- 
ing from a high-temperature initial state, the magneti- 
zation, its second moment and the autocorrelation have 
been measured at the critical point. Similarly, we have 
studied the behaviour of the magnetization, its deriva- 
tive with respect to the temperature and the Binder cu- 
mulant for a completely ordered initial state. Theoret- 
ically one expects a power-law for all these observables 
within the short-time regime at the critical point. This 
is indeed observed in our numerical simulations. All the 
dynamic exponents and static exponents as well as the 
critical temperature have been determined. The results 
are consistent and strongly support a full dynamic scal- 
ing in the short-time regime of the dynamic evolution. 
The values of the dynamic exponent z and the static 
exponents P/v and v are independent of initial condi- 
tions and agree well with those measured in equilibrium. 
Unlike non-local cluster algorithms, the short-time dy- 
namic approach studies the dynamic properties of the 
original local dynamics. Our measurements of the static 
exponents are comparable with large scale simulations in 
equilibrium. All this indicates that the dynamic mea- 
surements of the critical exponents are promising. The 
present work could possibly be extended to the diluted 
or the random field Ising models. 



Note added 

Recently, the results v = 0.6298(5) 
0.518(1) [H have been reported. 



and [3/v 
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